arXiv: 1505.03141v2 [hep-th] 31 Aug 2015 


Preprint typeset in JHEP style - HYPER VERSION 


TAUP-2994/15, FPA-15/10 


Parity Breaking Transport in Lifshitz Hydrodynamics 


Carlos Hoyos * 1 , Adiel Meyer 2 and Yaron Oz 2 

1 Department of Physics, Universidad de Oviedo, Avda. Calvo Sotelo 18, 33007, Oviedo, 
Spain 

2 Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, 
Tel Aviv 69978, Israel 


Abstract: We derive the constitutive relations of first order charged hydrodynamics for 
theories with Lifshitz scaling and broken parity in 2 +1 and 3 + 1 spacetime dimensions. In 
addition to the anomalous (in 3+1) or Hall (in 2+1) transport of relativistic hydrodynamics, 
there is an additional non-dissipative transport allowed by the absence of boost invariance. 
We analyze the non-relativistic limit and use a phenomenological model of a strange metal 
to argue that these effects can be measured in principle by using electromagnetic fields 
with non-zero gradients. 
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1. Introduction 

Experiments suggest that quantum criticality is underlying many of the exotic transport 
properties of ‘strange metals’ like high T c superconductors, heavy fermion, organic super¬ 
conductors, etc. In most cases very little is known about the quantum critical theory or 
theories that are responsible for the observed behavior, including something as basic as 
their symmetries. Any clue in this direction would then be extremely valuable. 

One of the simplest models of quantum critical behavior is the z = 2 Lifshitz theory 
introduced in [1], Theories with Lifshitz symmetry and dynamical exponent 2 ; exhibit a 
minimal set of symmetries including time and space translations, space rotations and the 
anisotropic scaling t —> A z t, x l —> Xx l , which may be anomalous for some observables. 
Starting with [2], strongly coupled theories with Lifshitz symmetry have been extensively 
studied using gravity duals with the goal of identifying universal properties that would 
help to understand better the nature of the quantum critical point. 

However, Lifshitz symmetry can in principle be extended to larger groups that include 
Galilean boosts and for z = 2 also non-relativistic conformal transformations 1 . For z = 1 it 

1 There are other possible extensions for discrete values of z < 1 e.g. [3]. 
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can be extended to the full relativistic conformal symmetry. One of the motivations for the 
present work is to identify effects that distinguish between systems with and without boost 
invariance. A first step in this direction is to develop the hydrodynamic description of a 
theory with Lifshitz symmetry at finite temperature. This is well motivated as an effective 
description of electrons in the case of strange metals [4]. The interaction rate between the 
electrons is much larger than the interaction rate with the ions in the lattice or impurities, 
so the electrons have time to reach local thermal equilibrium between scattering events 
with the other elements of the system that produce momentum loss. Note that this is not 
true in ordinary metals, that can be described as a weakly coupled gas of quasiparticles. 

The description of Lifshitz hydrodynamics was initiated in [5] for neutral fluids and 
extended for charged fluids and superfluids in [6]. In the present work we extend the 
charged hydrodynamics to include parity breaking effects. This is partly motivated by the 
fact that in many cases the strange metal behavior occurs in a magnetic material and/or 
in the presence of magnetic fields. More generally, parity breaking introduces many effects 
in hydrodynamics such as Hall transport in 2 + 1 dimensions or anomalous transport in 
3 + 1 dimensions. An interesting question is if new effects are possible when there is 
no boost invariance. We find that this is indeed the case in 3 + 1 dimensions, where 
there can be currents in the direction of magnetic fields or vorticity. Usually these kind 
of transports are associated to the presence of chiral anomalies and have been studied 
in multiple systems, including e.g. astrophysics [7,8] quark-gluon plasma [9], condensed 
matter [10, 11], cosmology [12] and relativistic hydrodynamics [13-16], 2 where they are 
commonly known as Chiral Magnetic and Chiral Vortical Effects. The novelty is that in 
Lifshitz theories it is allowed even if the currents are non-anomalous. In 2 + 1 dimensions 
there are additional terms that affect the Hall conductivity and Hall thermal transport. 

The transport described above is in the context of non-invariance under Lorentz boosts. 
In addition, we take a non-relativistic limit to derive the Lifshitz hydrodynamics in the 
context of non-invariance under Galilean boosts. We use it as a phenomenological descrip¬ 
tion of the strange metal as a Drude model with a strongly coupled electron fluid, in the 
spirit of [17] for graphene. We take into account the effects of scattering with the lattice 
or impurities by adding a drag term to the hydrodynamic equations. Interestingly, we find 
that in 3 + 1 dimensions a Chiral Magnetic term survives in the non-relativistic limit. In 
2 + 1 dimensions there is an anomalous Hall effect in inhomogeneous configurations, i.e. a 
Hall current in the absence of magnetic held. 

The paper is organized as follows: in § || we review Lifshitz hydrodynamics and de¬ 
rive the constitutive relations in parity breaking theories. We constrain possible terms to 
first order in derivatives using the local second law of thermodynamics. This is done for 
theories without Lorenz boost invariance, but otherwise relativistic. In § || we take a non- 
relativistic limit and derive the constitutive relations for theories without Galilean boost 
invariance. In § |^ we study some phenomenological consequences of broken boost invari¬ 
ance using the Drude model of non-relativistic hydrodynamics with a drag term. Finally 
in § |5] we summarize and discuss our results. We have collected several technical results in 

2 In this context, anomalous transport coefficients were first computed in holographic models [18-21]. 
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an Appendix |A|. 

2. Lifshitz hydrodynamics with parity breaking terms 

The lack of invariance under Lorentz boosts in theories with Lifshitz symmetry implies that 
the energy-momentum tensor is not necessarily symmetric. Using rotational invariance one 
can show that no asymmetric terms are expected in the hydrodynamic energy-momentum 
tensor at the ideal level, and we confirmed this explicitly by a calculation in free held 
theory [22]. Note that rotational invariance is unbroken in the rest frame of the fluid. 
For an observer that sees the fluid moving with a velocity u^ the corresponding symmetry 
condition is 

(T^ - T^)P m P v p = 0 , ( 2 . 1 ) 

where P^ = rj 111 ' + u^u 1 ' is the projector in the directions transverse to the velocity. 
The absence of invariance under Lorentz boosts on the other hand allows the following 
components to be asymmetric: 


(T^ - T^)u^P v p + 0. (2.2) 

In the rest frame of the fluid this becomes 

T 0i _ T i0 Q (2.3) 

But in a state with rotational invariance the expectation value of T° l vanishes, so at 
the ideal order the hydrodynamic energy-momentum is symmetric. When the velocity or 
other fluid variables are not constant there can be asymmetric terms that depend on the 
derivatives of the fluid variables. In previous works [5,6] we found the possible asymmetric 
terms to first viscous order in fluids and superfluids with unbroken parity symmetry. Here 
we will give a detailed presentation for charged fluids with parity breaking terms. 

We will work in the Landau frame T^ u u v = —eu^, where the energy-momentum tensor 
takes the form 3 

= evPu v + pP^ + 7 t^ u) + vrjf 1 + (u^n [ A ] + u u n [ ^ ] )u a . (2.4) 

Both 7 ta and ns contain terms that are at least first order in derivatives. The last symmetric 
term depending on nA ensures that the Landau condition will be satisfied when we add 
the antisymmetric part. In addition, the symmetric part should be transverse ng’u v = 0. 
The condition ( |2.lD constrains the anisymmetric part to be of the form 

Trjf' 1 = u^V A ] . (2.5) 

Where V A u v = 0. Then, the form of the energy-momentum tensor in the Landau frame 
becomes 

= euPv? + pP^ + n^ v) + u^V v A . (2.6) 

3 We follow the usual bracket convention for symmetrized A^B V ' > = |( A^B 1 ' + A V B^), and antisym¬ 
metrized - A"B^) tensors. 
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In addition to the energy-momentum tensor, there is a conserved global current whose 
constitutive relation is 

J* = qu /* + iS, (2.7) 

where we impose the condition = 0, so that the charge density is q = 

The name of the game is to find all possible first order viscous terms, constructed from 
derivatives of the velocity, temperature T, chemical potential n and external electromag¬ 
netic fields A t j. We consider the derivatives of the external fields to be of first order. There 
are many such terms that one can build, but we will constrain them by demanding a local 
form of the second law of thermodynamics holds. This means that it should be possible to 
define an entropy current (which coincides with the usual notion of entropy at equilibrium) 
whose divergence is not negative 

> 0 . ( 2 . 8 ) 

Starting with Landau, this approach has been used with great success to constrain hy¬ 
drodynamic coefficients. For parity breaking terms we will restrict the analysis to 2 + 1 
dimensions and 3+1 dimensions. If parity is unbroken the entropy current takes the generic 
form 

f s = ^ - ^V. (2.9) 

Where s is the canonical entropy density e + p = Ts + [iq. When parity is broken it is in 
principle possible to add to the entropy current additional terms that are pseudovectors, 
we will denote them as Va ■ The entropy current is then defined as 

# = su< 1 (V + D n Va ■ (2.10) 

We can derive the equation for the entropy current from the conservation equations 

d^ v u v + + i-iA. (2.11) 

Here we allow for an anomaly = A in the global current, such term can be present in 

3 + 1 dimensions. In 2 + 1 dimension one should simply set A = 0. Expanding explicitly 
the equation above one finds 

—Td^(su^) + d^TTg 1 ^ u u + u^d^V^u u + iid^ = -E + /iA. (2-12) 

Where we have used e + p = Ts + /iq, the first law 

dfj,£ = Td fl s + lid^q, (2-13) 

and we have defined the electric field as E^ = F^ u u v . Straightforward algebra leads to the 
following equation for the entropy current 

d,j': = rirV, - + d^DaVf). (2.14) 

Where a^ = u a d a u 11 is the acceleration and we have defined the vector 

= E tt - TP”d v (£) . (2.15) 
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The positivity condition d p js > 0 implies that terms on the right hand side of the equation 
that do not cancel out should form squares. We can separate parity preserving and parity 
breaking contributions 

dfijs = Ap + Ajj>. ( 2 . 16 ) 

They do not mix, so we should impose the positivity condition separately for each of them 


A P >0, > 0. (2.17) 

For illustration purposes we will review the analysis of [ 6 ] for the parity preserving con¬ 
tributions. The most general form of the viscous terms allowed by the entropy current 
equation are 

74 ^ = - rja^ - (P^d Q u a , 

V% = - axoP - 2a 2 £ fi , (2.18) 

v 11 =2a 3 a M + oE». 


Where r] and C are the shear and bulk viscosities respectively, a is the conductivity, and 
the shear tensor in d + 1 dimensions is defined as 


= P m P^ ( 2 d {a u p) - jP a0 (d a u^ 


(2.19) 


From ( [2.14 ), the positivity condition imposes the following constraints on the parity pre¬ 
serving coefficients 

77 > 0 , C>0, cr > 0 , aa\ > {a 2 + a^) 2 . ( 2 . 20 ) 


The combination a 2 — a 3 drops from the divergence, so it is dissipationless. It will vanish 
if we impose Onsager’s relations. 

We will now proceed with the parity breaking terms. The analysis depends on the 
number of dimensions, so we treat them separately. 


2.1 Parity breaking transport in 3 + 1 dimensions 

In 3 + 1 dimensions there can be interesting new effects when parity is broken due to 
chiral anomalies. A purely hydrodynamic analysis shows that the value of some transport 
coefficients is fixed by the anomaly [13-15] 

There are two possible independent pseudovectors we can construct [13], the vorticity 
and the magnetic field B They are defined as follows: 


'J l = \e^^u v d p u a , 

= l^t>° UvFpa . 


( 2 . 21 ) 


The entropy current has then the form 


% = su» - ^ + D^ + D b B>\ 


( 2 . 22 ) 


- 5 - 





The parity breaking contribution to the divergence is 

+ d,(D^) + d^DsB*) - £ CE, (2.23) 

Where we have introduced a possible chiral anomaly for the current A = CE^B^. Note 
that there are no u 2 or B 2 , so the positivity condition for the entropy current is that there 
should be an exact cancellation among the parity breaking contributions 


A^ = 0. 


(2.24) 


This implies that all the related transport is non-dissipative. The divergence of the vorticity 
and of the magnetic held are not independent scalar quantities, but they can be decomposed 
as follows: 

= 2a M uA = a, U B» - 2u^. (2.25) 

The possible parity breaking contributions that can appear in the viscous terms are 


V a^ = ~ T/W* - TP b B», 
^ 


(2.26) 


It is sufficient if the cancellation is realized on-shell, this means that at this order we are 
allowed to use the hydrodynamic equations truncated to ideal order. In particular, this 
fixes the acceleration 

n 1 

-P^duP. 


(A = 


q -E* - 1 


£ +P 


£ + p 


Collecting terms, we find 


A^ = ^ (d»D u - (2 D u + p u ) - Zudu ( 


+ (2£> w + ft,) 


+ B» d,D B - 


dyP 

£+p 


- 2 D b + ^ 
£ + p T 

(D b + /?b) - 


+ B P E,, 


C 1 {Db + / 3 b ) + Cb _ C JJ_ | = 0 

e + p T T 


(2.27) 


(2.28) 


Each of the brackets has to vanish independently, this results in the following conditions: 


d„D u - (2 D u + ft,) = 0, 

£ + p Vi / 

(2 D u + ft,) -2D b + ^ = 0, 

£ + p T 

d lt D R - ( D B + p B ) - Cb/3/4 = 0, 

e + p T 

q ( P B + /3 B ) Cb _ c P_ = „ 

£ + p T T 


(2.29) 

(2.30) 

(2.31) 

(2.32) 
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The derivation of the solutions can be found in the Appendix |A.1| , here we present the final 
result: 


PUfi,T) = 2jlTp B + b U] , (2.33) 

£b(/2,T) = C (» - - -^(d B + (5 B ), (2.34) 

V 2 £ + p J £+p 

U ,,T) = c(S- 2 -*f) + 2T(d B ~^Li B ) 

V 3 £+p) V e +p J 

- (2 d w + 2Tp/3 B + b w ), (2.35) 

£ + p 

D B (p, T ) = -CTp 2 + d B , (2.36) 

T) = -CT 2 p? + 2 JlTd B + d u , (2.37) 

d B (T) = 7 B T + T J T t^p-dr, (2.38) 

du (T) = 7 J 2 + T 2 J ^ppdr. (2.39) 


With b u (T) and (3 B {T) arbitrary functions of the temperature. 

For z = 1, we would expect (5 B = c B T and b u = c u T 2 if there is scale invariance, 
as in theories with Lifshitz symmetry. However scale invariance would actually be broken 
because for this temperature dependence when we do the integrals over temperature in 
(2.38), ( 2.39D we find logarithmic terms 


d B = 7 B T + 7 bT log(T/A), d u = 7 JT 2 + c u T 2 log(T/A), (2.40) 


where A is the scale associated to the breaking of symmetry. Therefore, Lifshitz symmetry 
demands (5 B = b u = 0. In this case the transport coefficients are simplified to 


/) = 0 , 

£ B (fi,T) = C ( p- 


1 


Qd 


qT 2 


2 £ + p) £ + p 


IB, 


2 2 Qd 


Up, t ) = c ^-- 3£+p 

D B {p,T)= l -CTp 2 + lB T , 


+ 2 1b T 2 - (2 1b[ iT 2 + 7 w T 3 ) , 

£ + p 


= -CT 2 p 3 + 2fiT 2r y B + ^T 2 . 


(2.41) 

(2.42) 

(2.43) 

(2.44) 

(2.45) 


Note that these are the same results as for the relativistic theory [18]. If we further impose 
CPT invariance (note that for a Lifshitz theory this is not a requirement), we will fix 
7 ^ = 0. For a detailed explanation on how the hydrodynamic quantities transform under 
CPT see [23]. 

2—z 2 

For z / 1, we would expect [3 B = c B T * and b u = c u T ~ if there is scale invariance, 
as in theories with Lifshitz symmetry. However, in this case the terms proportional to 75 
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and ju break scale invariance. In this case, Lifshitz symmetry demands 7 b = 7u> = 0. 
Therefore, the transport coefficients are 


/M/b T ) = (ZcbH + Ccj) T~* 


£,b(h,T) = C | n - 


1 qn 2 


Ufi,T) = Ci v- 


2 e + p 

3 


- Cb 


2 2 QR 


+ Cb : 


2 - z q 
2(1 - z) £ + p 

z 


T~*, 


3 £ + p) 1 — z 
+ 2 csji ) T" 


i_i« ir? 

£ +P 


£ + p \1 — Z 

Db^T)= 1 -CT[£ + -£—c b T^, 

Doj(fl,T) = -CT^fl 3 + _ ^ (2/Icb + Co,) T = 


(2.46) 

(2.47) 

(2.48) 

(2.49) 

(2.50) 


2.2 Parity breaking transport in 2 + 1 dimensions 

We follow closely the analysis of parity breaking fluids in [24], There are three possible 
independent transverse pseudovectors we can construct 


Ui =e^ p u v a p , 

U p =e*">u v E p , (2 . 51) 

U p = -^ p u u £ p . 

We can also construct two independent peudoscalars, the vorticity to and the magnetic field 
B, and a pseudotensor, the Hall viscosity a pu : 4 

u =^e' J ' vp u p d l ,u p , 

B =L^ p Ufi F vp , (2.52) 

=- A { e ^u p a: + e^u p a p ). 

Note that our definitions differ from those of [24], the relation is 7? t h e irs = — -B ours and 
^theirs = — 2cj ours . It will be convenient for the analysis of the entropy current to group all 
the possible pseudovectors (transverse and not transverse) in the following basis (R = j+^) 

V( l =e pup u u d p T = —TU P - RT 2 U P , 

V p =U P , 

vf =^ p n w d,£ = US + if 7S, (2.53) 

vt = US - Bu“, 

v p =e^ p d u u p = -U p - 2w/. 

4 An useful relation is 2d[ p u v ] = —2e llvp u p u}. 



We can add the following parity breaking terms in the viscous contributions to the 
energy-momentum tensor and the currents: 

• Symmetric part of the energy-momentum tensor: 

- CbP^B - (2.54) 

• Asymmetric part of the energy-momentum tensor: 

3 

= (2-55) 

i= 1 


• Current: 


• Entropy current 


i= 1 


j 1 : = su> i - 

i =1 


(2.56) 


(2.57) 


The parity breaking contribution to the divergence of the entropy current is 

5 


- ^t/W - + £ I 'gid.Tvf + ^B^vf + Wj? 


1 


i= 1 


di>i 


dv; 


dT 


3/x 


(2.58) 


Where we have set the anomalous term to zero .4 = 0 since the number of spacetime 
dimensions is odd. The positivity condition for the entropy current is that there should be 
an exact cancellation among the parity breaking contributions, since it is not possible to 
form squares 

Aq = 0. (2.59) 

This implies that all the related transport is non-dissipative. The divergences of the pseu¬ 
dovectors are not independent scalar quantities, but they can be decomposed as follows: 


= -2 um Q d a T - Ufd^T, 

(2.60) 

d/jV^ = d a (Bu a ), 

(2.61) 

= -2 um a d Q ^ - Ufd^, 

(2.62) 

& 

II 

o 

(2.63) 

= 0 . 

(2.64) 


It is sufficient if the cancellation is realized on-shell, this means that at this order we are 
allowed to use the hydrodynamic equations truncated to ideal order. In particular, this 
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fixes the acceleration to (2.27) and we have the following relations among derivatives 


=-!(!)/• 

P^ a d Q T = - TU? - RT 2 Ug = -Ta^ + RTS^, 
P^daU =U>i + = i(-^ + Efa, 

P^ a d a p = - (e + p)(a 11 - RE**). 


(2.65) 


Collecting terms, we find 


Ap = U 1 U 3 (fa - RT 2 fa - T5 x - RT 2 {d T fa + fa) + (dpfa + fa) + T(dpi) 1 - d T fa)) 
+ U\U 2 (fa + ^fa + dflUO r^ - d T fa - Tdrv^j 

+ U 2 U 3 (fa + (pfa + RT(d T i >3 - dpfa) - dpfa + RT 2 d T v^\ 


+ B6 


1 


-Xb + T 



drfa + 
q 



+ 2^6 I —fa 


+T ( I ), (ar " 5 ++ f (|), ^ + " 3) 


+ d^faU-2 + faBO + z fau^d^B 


( 2 . 66 ) 


Each of the brackets has to vanish independently, 


this results in the following conditions: 


z> 2 = 0 (2.67) 

fa - RT 2 fa -T 81 - RT 2 (d T fa + fa) + {dpfa + fa) + T(dp,fa - d T fa) = 0 (2.68) 

fa + j;fa + + dp,fa - d T fa - Td T fa = 0 (2.69) 

fa + j;fa + RT{d T fa - dp fa) - dp fa + RT 2 d T fa = 0 (2.70) 

^ Cw+t (^) ( 5 tz> 5 +z>i)+ KS) ( ^ 5+z>3)= ° (2 - 72) 


For fa = 0 we recover the results of [24]. If we assume Lifshitz scale invariance, in general 
the coefficientes will have a power dependence on the temperature times an arbitrary 
function of / 2 . The power is determined by demanding that all the terms in the same 
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equation have the same scaling, taking into account that the scaling dimension of the 
temperature is [T] = 2 . This fixes: 


1 — 2 
T—, 

+2 

1 — 2 

^ z 

i 

1 — 2 


1 — 2 

~ 1 

T—, 

-52 

~ T z , 

$3 ~ Tz, 

1 — 2 
T~, 

h 

~ T* , 

1-2 _ 1 
Z /4 ~ T 2 , Z /5 T z 


C B ~ T^ Z , Co; ~ T z 


(2.73) 


3. The non-relativistic limit 


We now study fluids with broken Galilean boost invariance. Instead of deriving the hy¬ 
drodynamic constitutive relations from the start, we will take a non-relativistic limit of 
the Lifshitz hydrodynamic equations we have discussed in the previous sections. This will 
automatically lead to non-relativistic hydrodynamic equations. Similar limits have been 
used in the context of Galilean invariant fluids in [25,26]. 

In the relativistic fluid the maximal velocity c appears in uJ l = (1, /3) /yA — /3 2 , where 
/+ = v l /c. In the non-relativistic limit c —>• oo, the pressure is not affected while the 
relativistic energy is expanded in terms of the mass density p and the internal energy U as 

e = c 2 p-^ + U (3.1) 


From the thermodynamic relation e + p = Ts + pq, we get that 

v 2 Pnr 

q = pc -p p = c -\ - 


(3.2) 


The electromagnetic fields scale as A % —> Aj, Aq —>• A t /c, which give the relations: 
In 3+1 dimensions: 


E l — — (^E l NR + e l -i k VjBN Rt k'j + o(l/c), B l — B l NR + 0(l/c 2 ). 


(3.3) 


In 2+1 dimensions: 


E l = - c {E l NR - & Vj B NR ) + o(l/c), B = B nr + 0(l/c 2 ). (3.4) 

The subscript NR will be left out in the reminder of the text. We use the conventions in 
which e 012 ' 3 = e 123 = e 12 = 1 . 

For charged hydrodynamics we first go to the Eckart frame, where = qu + To first 
order in derivatives the energy-momentum tensor in this frame is 

= evAvA + pP^ + u»VX + tA 1 ' - + u u ^). (3.5) 


Note that in the Eckart frame the current conservation equation gives the usual continuity 
equation. If v l ~ 0(1), the energy conservation equation becomes 


d t p + di(pv l - v l ) = 0 . 


(3.6) 
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This would make the charge and mass currents different, while we are interested in de¬ 
scribing systems made of particles with fixed charge over mass ratio. In order to have a 
consistent expansion, the dissipative terms in u l should start contributing only to the next 
order, where they enter in the thermal current. Then we have to impose that v i 0(1/ c 2 ). 
In this case, there are no contributions from v l terms to the Navier-Stokes equations. 

Galilean boost invariance will fix a relation between momentum density and mass 
current P l = J l = pv l . When we take the non-relativistic limit of the boost-breaking 
hydrodynamics, this condition does not hold anymore beyond ideal order. We find that 
the Navier-Stokes equations become 

• In 3 + 1 dimensions: 

d t P i + dj (PV) + &p = 
p (/;' + ( ijk v 3 Ih) + dj (ria ij + &Sd k v k ) . 

Where cr^- = diVj + djVi — ijdkV k is the shear tensor (in this case for d = 3). The 
momentum density is 

P l = pv i - a 0 a l - oit&T - T$ b B\ (3.8) 

where we define the acceleration as a 1 = Dtv 1 = dtv 1 + v 3 djv l and the magnetic held 
as B l = i e^ k Fjk. The term f3 B allows a Chiral Magnetic Effect in a non-relativistic 
theory. If the theory was obtained from a non-relativistic limit, it is determined by 
the value of the anomaly in the original relativistic theory. Note that the continuity 
equation for the current holds in this limit, and the effect of the anomaly enters in 
the thermal current and the momentum density. 

We got these expressions by demanding that the leading order dependence of the 
relativistic coefficients with the speed of light is, 

r/~ 0 (c), C ~0(c), oi ~ O (c 3 ) , /3u, ~ O (c 2 ) , /3 B ~ O (c), a 2 ~ O (1) , 

a 3 ~0(l), a ~ O ~ O , i B ~ O , C ~ O . (3.9) 

In principle we could add to the momentum a term of the form 


PJti, (3-10) 

where II* = is the vorticity. However, the entropy analysis that we detail 

in Appendix [A.‘2| reveals that those terms should be zero in the non-relativistic limit. 
Therefore, Chiral Vortical Effects are suppressed in the non-relativistic limit. 


• In 2 + 1 dimensions: 

d t P L + dj (PV) + Wp = 

p (. E 1 - e lj VjB ) + dj + r]Hd lj + (dk,v k ^) 


12 


Where the Hall viscosity tensor is 


a'i = (3.12) 

and &ij is the shear tensor for d = 2. The momentum density is 5 * 

P ! = pv l - Uatf - a T d i T - (i T ^d,fr - ^d jf Jt NR - P E j j (Ej ~ e jk v k B). (3.13) 
We got this expression by demanding that the leading order dependence of the rela¬ 


tivistic coefficients with the speed of light is, 

(u ~ Vh ~ 0(c), Cb ~ 0(1), (3-14) 

h = |#3 0) + c Ai 1) +0(c), M3 = c 2 /4 0) + P? + 0 ( 1 ), M2~0(c 2 ), (3.15) 

5i = 5[ 0) + 0(l/c 2 ), 5 3 = + 0(l/c 3 ), 82 ~ 0(l/c). (3.16) 

In principle we could add to the stress tensor terms of the form 

6 ij (( B B + ( n n), (3.17) 


with B the magnetic held and = —^e^diVj the vorticity. However, the entropy 
analysis that we detail in Appendix A^ reveals that those terms should be zero in 
the non-relativistic limit. 


4. Drude Model 

We model the collective motion of electrons in the strange metal as a charged fluid moving 
through a static medium, that produces a drag on the fluid. The hydrodynamic equations 
are 


d,,,! 11 = A, = F va J a -\c5 ui J i . (4.1) 

We are interested in describing a steady state where the fluid has been accelerated by the 
electric held, increasing the current until the drag force is large enough to compensate for 
it. In order for this to happen we keep the external helds constant in time, although we 
allow them to vary slowly in space. In the steady state conhguration the current does not 
change but other quantities like the energy can change with time since there is dissipation. 

4.1 Effects in 3+1 dimensions 

The Navier-Stokes equations with the drag term included are 

d t P l + dj (PV) + d V = 

p (e 1 + + dj ( w ij + S ij Cd k v k ) - A pv\ 

5 The coefficients /? are combinations of the original coefficients fl in the non-relativistic limit 

separated the different independent contributions. 


(4.2) 

. We have 
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Where P l is given by (|3.8| ) . We can see the effect of breaking of boost invariance and parity 
by introducing constant electric and a slowly varying magnetic field E x and B z (x), in such 
a way that when we solve for the velocity derivative terms are suppressed. We also assume 
that the amplitude of the electromagnetic field is small, so non-linear terms are subleading. 
To leading order, the current has only a longitudinal component 

J x = —Tlx- (4-3) 

Since d{V 1 = 0 to this order, we can take the density p to be constant, as well as the 
temperature and the chemical potential. 

In addition to subleading corrections to the longitudinal current, currents in the trans¬ 
verse directions are also generated. There are two kind of contributions, pointing in different 
directions. One is coming from the Lorenz force term 

J y = ( 4 - 4 ) 

The second is due to the Chiral Magnetic term and points in the direction of the magnetic 
field 

TBr 

Jz = -^-E x d x B z . (4.5) 

Note that this current would be forbidden in a Galilean-invariant theory (with a single 
species of particles). 


4.2 Effects in 2+1 dimensions 

Assuming a time-independent configuration, the hydrodynamic equations for the density 
and the velocity are 


di{pv l ) = 0, d k {P l v K ) + d'p = d k [ya Kl + find* 1 + (5 kl e) - A pv l + p[E l - e' 3 Vj B). (4.6) 
Where the momentum P®, given by (|3 .13), is to first order 


P l = pv l - a a v k d k P - a T d l T - fa&djT - B^ 3 d jt i NR - Be^ (Ej - e^v k B). (4.7) 


We will allow a constant electric field E x and gradients of temperature and chemical po¬ 
tential along the same direction T(x), p(x). This implies that the pressure and the density 
are not constant but also vary along the x direction p(x), p(x). We will assume that the 
gradients are small so we can do a derivative expansion and solve order by order around 
constant density and velocities. 

To leading order we get the usual longitudinal current 

Jx = jE x . (4.8) 

A transverse (Hall) current is also generated. The leading contribution to the velocity in 
the y direction is 

-Be 


1 


|2 


|2 


|2 








This can be interpreted as an anomalous Hall effect: there is a transverse current in the 
absence of magnetic fields. 

There can also be a higher order correction due to the Hall viscosity O(c^). Therefore, 
this non-homogeneous configuration allows to measure the new transport coefficients from 
a Hall current generated in the absence of magnetic fields. In the case of a fluid with 
Galilean invariance the Hall viscosity would be the leading order contribution. One can 
distinguish them in principle because the Lifshitz contributions would be largest close to 
a minimum or a maximum of the pressure, temperature and or chemical potential, while 
the contribution from the Hall viscosity would approximately vanish. 


5. Summary and conclusions 


We have studied the constitutive relations of fluids for systems with Lifshitz symmetry and 
broken parity. We find new possible effects in the presence of magnetic field or vorticity. 
These are qualitatively different from theories with unbroken boost invariance and could 
be used to discern whether quantum critical points are really Lifshitz or a different type of 
theories with scaling symmetry. 

When the condition of boost invariance is relaxed there can be new terms in the 
energy-momentum tensor that can be grouped in a vector V£ 

T^ =eu?u v + pP^ + 7r { <r } + tfVX, 

5 (5.1) 

A =qvP + zA 


In 3 + 1 dimensions the terms that break parity are proportional to the magnetic field or 
the vorticity 


V A,p = - 77W - TPbB*, 

^ 


(5.2) 


If there was boost invariance, both and will be present only if there are chiral 
anomalies. In a Lifshitz theory this is not necessary, however in order to preserve the 
scaling symmetry for z = 1, one should set the coefficients (3^ = /3b = 0 . 

In 2 + 1 dimensions there are many more terms allowed 


A’ = - " C bP^B - w#*, 

v w = ~ T ilhv?, 

i—1 

i— 1 


Where the definition of the pseudovectors Vj 1 can be found in ( 12.51 ) , ( 2.53 ). They are 
constructed with epsilon tensors and the gradients of temperature and chemical potential 
and the acceleration and electric held. The entropy analysis imposes some constraints 
among the coefficients but all the terms are allowed. 
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In the non-relativistic limit the breaking of boost invariance can be seen as a modifica¬ 
tion of the momentum density, which is no longer equal to the current J l = pv l but receives 
corrections to first order. In 3 + 1 dimensions there is a parity breaking term proportional 
to the magnetic field 

P i = pv i - a a a l - a T d l T - TfoB*. (5.4) 

In this case the second law forbids a term proportional to the vorticity but the scaling 
symmetry does not impose the condition that /3b = 0. Therefore the relativistic and non- 
relativistic cases are qualitatively different. Using a Drude model with drag coefficient A 
for the strange metal, the parity breaking term is responsible for producing a current in 
the direction of the magnetic field, if the field is inhomogeneous 

d x B z E x . (5.5) 

In 2 + 1 dimensions the momentum density receives several parity breaking contribu¬ 
tions 

P' = pv l - a a v k d k v l - ot T d l T - /3 T e lj djT - d 3 p, NTt - /3 E e lj {Ej - e k v k B). (5.6) 

Using the Drude model, gradients of the temperature and/or chemical potential will intro¬ 
duce an anomalous Hall conductivity 

J y = i —dip ~ Pt 9 2 x T - fin dlp NR E x . (5.7) 

A l P 
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A. Entropy calculations 

A.l Positivity conditions in 3 + 1 dimensions 

In order to solve the conditions 

d»D u - (2 D u +p u )*teL- (£) = 0 , 

e + p \T J 

(2 D u + p u ) -2D b + ^ = 0, 

d^B ~ (D B + /3b) - ? B dA = 0, 

£ + p T 

q(B>B + /3b) Cb _ = „ 

e + p T T 


(A.l) 

(A.2) 

(A.3) 

(A.4) 
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We first take p = p/T and p as independent thermodynamic variables, so that 


dfiP B ,u — 


/ 9Db,i 


dpP + 


( dP B ,u 

V dp 


d+P- 


Collecting terms of the same type we get 

dD u \ 2ZA, + f3 u 


dp 

dDu 

dp 


£ + p 


= 0 , 


£,uj — 0 , 


(2 P u + 0 U ) -S— - 2D b + % = 0, 
E+P T 

( dP B \ _ P B + Pb _ „ 

\ dp ) - £ + p 


dP E 

dp 


- £b = 0 , 


c iP]t±M + tB_ c o. 
E + p T 


We can trade derivatives with respect to p by derivatives with respect to T by 
equations 

'dT\ qT 2 

dp ) P ~ 


dT\ 
dp). 


T 

E+P ’ 


E+P 


(A.5) 

(A.6) 

(A.T) 

(A.8) 

(A.9) 

(A.10) 

(A.ll) 
using the 

(A.12) 


Then, 


dP u 

dT 


dP, 


( 2 Du + Puj) 

++), 


E+p 

Combining the first three and the last three equations we get 

dl)j_, 


— (2P U + j3cj) — 0, 

(A.13) 


(A.14) 

~^r - idb + % = o, 

E + p T 

(A.15) 

— ( Db + Pb) = 0, 

(A.16) 


(A.17) 

] - + ^L-Cp = 0. 

(A.18) 


dp 


- 2 TP b = 0, 


( Wr- CTP =°■ 


(A.19) 
(A.20) 
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The last equation can be integrated trivially 


DsfaT) = i CTp 2 + d B (T ). 

Then, from the other equation we get 

DjjiX) = ^CT 2 fL 3 + 2fiTd B (T) + d UJ {T). 
From (A.1C) we get the equation 

T d' B — d B = P B . 


(A.21) 


(A.22) 


(A.23) 


This 


is is consistent only if = 0, hence (3 B only depends on the temperature T. 


Then 


v / 7 

the equation can be integrated and the general solution is 


d B (T) = 1b T + T J T ^ 

From ( A.13[ ) we get the equation 

Td' u - 2d w + 2jdT/3 B = f3 u 

This is a consistent equation only if 

/% = 2pT/3 B + b u (T). 

The general solution for d u is 

rT b u (r) 


dr. 


du{T) = lui T 2 + T 2 J ^-dr. 


Finally, we get the following expressions for £ B and £ u 

qT 


c n I 1 qfi2 
= C \ n - - 


2 £ +PJ £ + P 


(d B + Pb) 


(A.24) 


(A.25) 


(A.26) 


(A.27) 


(A.28) 


2 — 3 


= C ( p 2 - — ) +2 T[d B - ^-d B ) - -^-(2 d„ + 2Tfip B + b w ). (A.29) 

V 3 £ + pJ \ £ + p J £+p 

A.2 Non-relativistic limit in 3 + 1 dimensions 

The relativistic hydrodynamic equation ( |2.11| ) reduces to the non-relativistic form 

d t p + di (pv l ) = 0 (A.30) 

d t U + di (Uv l ) - p NR ( d t p + di ( pv *)) +pdiV l = rja 2 + C (dy) 2 + a [a 1 ) 2 

+ Ten + f3 B B l + ^a 2 &T S j + di (2a 3 a 4 + £ d l T + + &B*) - CE l Bi. 

(A.31) 
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We have defined the vorticity as 12* = ^e^ k djV k - In this expression we have already 
extracted the leading order dependence on c of the relativistic coefficients: 

V ~ O (c), C ~ O (c), a ~ O (c 3 ) , Aj ~ O (c 2 ) , Pb ~ O (c ), o 2 ~ O (1), 

a 3 ~0(l), cr~0^—V ^ ~ O (-1 , C ~ O (-V (A.32) 


In order to find the non-relativistic constraint equations which guarantee positive definite 
entropy current, we take the non-relativistic limit of (2.2S|), 


A p = Q 


+ 0 


12* ( diD w - (2+ /3 U ) — - UX 


£ + p 


P 


q 


+ O (1) B l ( diD B - 

'1 


12 l Ei (2+ 0 U ) —*— - 2Z2 S + 

e + p T 

dip 


(Db + Pb) - £,b9i 
e + p T 


+ 0 


B l Ei 


g{PB + Pb) + £b _ C P_ | = Q 
£ + p T T, 


(A.33) 


Using the first law of thermodynamics dU = Tds + pNRdp and U + p = Ts + p-nrP, 
Eq. (A.31) determines the divergence of the canonical entropy current in the non-relativistic 
case. This equation was derived by taking the 0(l/c) order of ( [2.1 ip , thus we demand 
that Eq ( A.33| ) will be of the same order. This fixes the leading order dependence of 
to be of order O (c°) and of Db to be O (-). After factoring out the c dependence and 
taking the limit, we find the following conditions: 


diD u - —d iP + i^diT = 0 
p T z 

(A.34) 

y = 0 

(A.35) 

d i D B -—d i p + ^d i T = 0 
p T z 

(A.36) 

Pb~% = 0 

(A.37) 


A.3 Non-relativistic limit in 2 + 1 dimensions 


The relativistic hydrodynamic equation ( 2.11 ) reduces to the non-relativistic form 

d t p + di (pv l ) = 0 (A.38) 

d t U + di (Uv l ) - pNRdtP ~ PNRdi (pv l ) + pdy 
C 


= W 2 + ^ {dy) + a (a*)" + ( w d k v k 12 + ( B Bd k v k + - (£3 - T 2 pi) Vjf 4 


+ Tfo&ViLj + ryyvy 


Pnr 


+ di[ - Si \T i + 5 2 e ij L 1 + S^d, 


Pnr 


(A.39) 


Where we have defined 


T* = - e ij djT , V = E l - e ij Vj B. 


(A.40) 
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In order to find the non-relativistic constraint equations which guarantee positive def¬ 
inite entropy current, we take the non-relativistic limit of (2.61). The result is 


A n = 0 (As - RT 2 ih - TS i - RT 2 (d T h + v i) + (d p h + u 3 ) + T{8 -^i - 5 r z> 3 )) 

+ 0(^-1 ( M2 + ^3 + ^^ 5 t + Ui + dgifi - drh - Td T i>4 


+ O ( 1 ) ( d 2 + + RT(d T i>3 - 9,-i^i) - dgiUi + RT 2 d T Vi 


+ o 

+ o 


1 


;Cb + T 


dp 

Y 


<9jT 4 + 


1 f dp 


T \dq 


2 T 


Cw + T 


<9p 

57 


(5t^ 5 + f'l) + 




1 f dp 


T \dq 


{dp.v 5 + £ 3 ) 


(A.41) 


Using the first law of thermodynamics dU = Tds + PNR^q and t/ + p = Ts + PnrQ, 
Eq. (A. 31) determines the divergence of the canonical entropy current in the non-relativistic 
case. This equation is of order 0(l/c), thus we demand that Eq ( A.41 ) will be of the same 
order. 


Given the expansion of the relativistic chemical potential with c (3.2), in the non- 
relativistic limit we can replace thermodynamic derivatives to 


dp, -* — d mR (A.42) 

8t —t d T + — (l + ^) dfj. N R (A.43) 

Therefore, using also that R = 4 + 0(c -3 ), 

- {u + p) 4 0) + /4 1} - r 2 AS 1} - Thi 0) - T 2 5 t^ 5 - T 2 z>i + h + T 2 d mR h, 

- Td T v 3 - d mR v 3 = 0, (A.44) 

+ ^3 0) = °> ( A - 45 ) 

d 2 + + Td T v 3 + 5 MjVi? i>3 - T ‘ 2 d llNR u\ + T 2 5 T i>4 = 0 , (A.46) 

0? = 0, (A.47) 

C, = 0. (A.48) 


The leading order dependence of the coefficients v with c is, 

~ O(c), v 3 ~ 0(1), z /4 ~ 0(1), P 5 ~ O(c). (A.49) 
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